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Abstract 

As the closest radio galaxy, Centaurus A is a pow- 
erful laboratory for the X-ray study of radio-emitting 
structures and their interactions with the hot interstel- 
lar medium (ISM). This paper details our interpreta- 
tion of the remarkable X-ray enhancement which caps 
the inner southwest radio lobe, at a radius of about 
6 kpc from the galaxy center. The shell of X-ray- 
emitting gas is hotter than the ambient ISM, and over- 
pressured by a factor of 100. We argue that it is 
heated compressed material behind the supersonically- 
advancing bow shock of the radio lobe, the first ex- 
ample of the phenomenon to be clearly detected. The 
results demonstrate that Cen A is actively re-heating 
nearby X-ray-emitting gas. The shell's kinetic energy 
is ~ 5 times its thermal energy, and exceeds the ther- 
mal energy of the ISM within 15 kpc of the center of 
the galaxy. As the shell dissipates it will have a major 
effect on Cen As ISM, providing distributed heating. 

1 Introduction 

There is much current interest in the possibility that 
radio sources heat the interstellar and intercluster 
medium. Such heating would help to explain the weak- 
ness or absence of lines from gas cooling below 1 keV 
in th e densest central regio ns of galaxies and clusters 
(e.g. JPeterson et all 12001). A good place for heating 
would be behind the bow shock of a supersonically ex- 
panding radio structure. In this paper we outline the 



simple theory and its application to Cen A, the first 
source to show clear and direct X-ray evidence of su- 
personic expansion and heating. 

2 Theory 

The standard model for the expansion of a power- 
ful radio source powered by a jet which is supersonic 
with respect to the X-ray-emitting interstellar medium 
(ISM) is illustrated in Fig. [J The jet terminates at the 
beam head (in a feature identified as the radio hotspot) 
where the jet fluid passes through a strong shock to in- 
flate a cocoon of radio-emitting plasma. The energy 
and momentum flux in the flow is normally expected 
to be sufficient to drive a bow shock into the ambient 
medium ahead of the jet termination shock. In the rest 
frame of the bow shock, ambient gas is heated as it 
crosses the shock to fill a region surrounding the lobe 
of radio-emitting plasma. 

The sound speed in gas of temperature T is given by 



IjkT 



where 7 is the ratio of specific heats, which for a 
monatomic non-relativistic gas is 5/3, k is the Boltz- 
mann constant, toh is the mass of the hydrogen atom, 
and /imn is the mass per particle, which is 0.6 run 
for a gas with normal cosmic abundances. Under 
these conditions, c s m 516(fcT/keV) 1 / 2 km s _1 rj 
0.54(feT/keV) 1 / 2 kpc Myr" 1 . 



1 
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The Mach number of the speed of advance, t> a d v > of the 
bow shock into the ambient medium is M. = f a dv/c s > 
which in convenient units can be expressed as 



core (nozzle) 




Ambient Gas 

Figure 1 : In the standard model for powerful radio sources, 
a supersonic jet (yellow) terminates at the beam head, pro- 
ducing a radio hotspot. Provided the shocked radio-emitting 
fluid forming the radio lobe (pink) has enough internal en- 
ergy or momentum density to drive a leading bow shock, 
ambient X-ray-emitting gas will be heated as it crosses the 
shock to fill the blue region. 



M « 580(v adv /c)(£;T/keV) 



-1/2 



(2) 



where c is the speed of light. In a simple application 
of the Rankine-Hug oniot conditions for a strong shock 
(e.g., Spitzei , 1978h . the pressure, density, and temper- 



ature, respectively, in the unshocked (subscript 1) and 
shocked (subscript 2) regions at the head of the bow 
shock are related by 



P 2 /Pi = (5M 2 - l)/4 



p 2 / P i = 4M 2 /(M 2 + 3) 



(3) 



(4) 



T 2 /Ti = (5M 2 - 1) (M 2 + 3) /16M 2 (5) 
for a monatomic gas. 

If the density in an X-ray-emitting gas is described by 
the proton number density, t] p , then the thermal pres- 
sure is given by 



P = rjpkT/Xn 



(6) 



where X is the abundance of hydrogen by mass, 
which is 0.74 for normal c osmic abundances (e.g., 
Birkinshaw & Worrall[[T993L . 



For hot X-ray emitting gas where the cooling in line 
radiation is unimportant, the X-ray emissivity, £ , be- 
tween energies E\ and E 2 depends on temperature and 
proton density approximately as 



2th0.5 



£ oc r)pT' 



-Et/kT 



-E 2 /kT^ 



(7) 



where the weak energy dependence of the Gaunt fac- 
tor is ignored. At temperatures below ~ 1 keV, line 
radiation cannot be ignored, and the plasma models 
incorporated into X-ray-spectral-fitting programs such 
as XSPEC can be used to find the dependence of the 
emissivity on energy. We find that in the energy band 
0.8-2 keV, where Chandra and XMM-Newton are 
most sensitive, for an A4 = 4 shock, the X-ray emis- 
sivity contrast between shocked and unshocked gas is 
a factor of 3 higher if the ambient gas is at 0.29 keV (as 
found for Cen A, Table Q than if the external medium 
has a typical cluster temperature of 4 keV. 
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Complications apply in reality, and in practice these 
are difficult to treat even with data from observa- 
tories as powerful as Chandra and XMM-Newton. 
Firstly, there is observational evidence that in super- 
nova remnants with shocks of comparable Mach num- 
ber to that found for Cen A the post-shock el e ctrons 
are cooler than the ions (e.g.. iHwangetall 2002; 



iRakowski. Gh avamian & Hughes], 2003b . This is not 
taken into into account in our modelling. Secondly, 
the simple Rankine-Hugoniot equations that we quote 
and apply do not take into account the fact that the bow 
shock around a lobe is oblique away from its head, with 
a consequent change in the jump condi tions and the 
emissivity contrast (e.g.. I WilliamsL 1991 ). However, if 
Cen A's shell represents a spherical expansion rather 
than a lobe structure, then the shock should be normal 
everywhere and the equations above will hold. 

3 Historical Perspective 

It is possible to interpret ROSAT-detected X-ray cavi- 
ties coincident with the inner parts of the radio lobes 
of Cygnus A as due to an emissivity contrast be- 
tween bow-shock heated gas outside the lobes (heated 
to temperatures above the ROSAT energy band) and 
the more easily detected am bient cluster medium 
( Carilli. Perlev & HarrisLll994b . However, the parame- 
ters of the shock are not effectively constrained by the 
ROSAT X-ray data. Similarly, Chandra observations 
of Cygnus A find gas at the sides of the lobes to have 
kT ~ 6 keV, slightly hotter than the value of 5 keV 
from ambient medium at the same cluster radius, pos- 
sibly indicating cooling after bow-shock heating, but 
again the data do not us efully constrain model parame- 
ters (Smi th et all l2002h . Other reports of lifting of gas 
(leading to eventual heating) by radio lobes or hot bub- 
bles are not thought t o involve supersonic expansion 
(e.g., Churazov et al. , 2001 : Ouilis. Bower & Balogh, 
200lb . although it has been suggested recently that fila- 
ments of hot gas in the atmosphere of M 87 are slowing 
from s upersonic speeds after ejection from th e galaxy 
center dForman et al.ll2004l:Eraft et alll2004b . 

The first and best example of a shell of heated gas 
which can reasonably be attributed to supersonic ex- 
pansion is in Cen A. High-quality Ch andra and XMM- 
Newton data JKraft et alll2003ll2004 provide the tem- 
perature and density constraints needed to test the 
model and measure the supersonic advance speed of 
the bow shock responsible for the heating. 



4 Cen A 

Cen A is our nearest radio galaxy, at a distance of 
3.4 Mpc ( Israeli 1998b so that 1 arcsec corresponds 
to ~ 17 pc, and is an example of a low-power radio 
galaxy. In such sources the radio jet is normally ex- 
pected to have slowed co nsiderably throug h entrain- 
ment of ambient material dBickneli fl994h . at which 
point the model described by Fig. ^ no longer holds. 
The full extent of Cen As radio emis s ion co vers sev- 
eral degrees on the sky Jjunkes et all Il993b . Within 
this l ies a sub-galaxy-sized double-lobed inner struc- 
ture (burns. F eigelson & Schreier . 1983b with a pre- 



dominantly one-sided jet to the northea st and a wea k 
counter-jet to the southwest (IHardcastle et all |2003), 
embedded in a radio lobe with pressure 1.4 x 10~ 12 Pa 
or more, greater than the pressure in the ambient ISM 
(~ 1.8 x 10~ 13 Pa; Table QJ, and so which should be 
surrounded by a shock. Around this southwest lobe 
there is a shell of X-ray emitting gas which appears to 
have the geometry of th e shocked ambient gas in Fig. [J 
( Kraft et all I2003L 120041) . Although the capped lobe is 
around the weak counterjet, so it is not evident that the 
lobe is being thrust forward supersonically with respect 
to the external interstellar medium (ISM) by the mo- 
mentum flux of an active jet, the current high internal 
pressure in the radio lobe ensures its strong expansion. 

5 Application of the model to Cen A 

The temperature, proton density and pressure of 
the ambien t ISM and the X-ray shell, taken from 
Kraft et all (120031) are given in Table The ambient 
medium is measured to have rj p ~ 1.7 x 10 3 m~ 3 and 
kT = 0.29 keV, whereas the shell is ten times hot- 
ter, at kT = 2.9 keV, and twelve times denser, with 
T] p ~ 2 x 10 4 m~ 3 . From equations (4) and (5) above, 
we see that temperature and density measurements for 
both the ambient medium and the shocked gas directly 
test shock heating, since only two of the four parame- 
ters are required to measure the Mach number, and the 
other two test the model. 

The most straightforward application of the equations 
finds an inconsistency, since the densities and temper- 
atures are not self-consistent. The shell's density and 
temperature are wrong for gas directly in contact with 
the bow shock. However, we can find a Mach number 
consistent with shocking the gas to a temperature and 
density such that the combined thermal and ram pres- 
sure is in pressure equilibrium with the thermal pres- 
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Table 1 : Physical parameters of the gas in various regions of Cen A 



Structure 


kT 


Vp 


Pressure 




0.4-2 keV 




(keV) 


Proton density 


(Pa)t 




relative 






Cm" 3 ) 






emissivity, £ 


ISM (measured) 


0.29 


1700 


1.8 x 10- 


13 t (thermal) 


1 


Behind bow shock (inferred) 


6.8 


6530 


2.1 x 10" 


11 (thermal+ram *) 


13 


Shell (measured) 


2.9 


20000 


2.1 x 10" 


11 (thermal) 


127 



f 1 Pascal =10 dyn cm 2 J incorrectly reported in Table 5 of Kraft et al. (2003) * Pi v 'tdv/^ 
By comparison, the minimum-energy pressure in the radio lobe is ~ 1.4 x 10~ 12 Pa. 



sure of the detected shell: A4 = 8.5, t> a d v ~ 2400 km 
s _1 . The post-shock temperature is kT2 ~ 6.8 keV. 
The 6.8 keV gas flows back from the shock, into the X- 
ray-detected shell at 2.9 keV. The characteristics of this 
undetected hotter gas are given in Tabled In this table 
we also quote estimates of the relative X-ray emissiv- 
ity (per unit volume) of gas in the different structures 
over the 0.4-2 keV energy band, where the Chandra 
response peaks and is relatively flat. The gas directly 
behind the bow shock has a predicted emissivity that 
is an order of magnitude fainter than that in the shell, 
accounting for its absence in our measurements. 

6 Conclusions 

• A hot shell of X-ray-emitting gas caps the south- 
west radio lobe of Cen A. 

• Shock jump conditions for an advance speed of 
~ 0.008c (M = 8.5) are satisfied if the shell 
is in pressure balance with unseen gas at kT ~ 
6.8 keV behind the bow shock. The emissivity 
of the 6.8 keV gas is too low to separate its X- 
rays from those of the ten-times-brighter gas of 
the shell. 



The radiative timescale for material in the shell 
(~ 2 x 10 9 yrs, using equation 5.23 of SarazirJ 
( 1986)) is large compared with the lobe expansion 
time (< 2.4 x 10 6 years), so the ma terial in the 
shell behaves as an adiabatic gas (e.g. JAlexa nder. 
2002). 



The cooling has improved the shell-ISM contrast 
by an order of magnitude, an important factor in 
the detection and modelling of the shell. 

The lobe is (or was) powered by energy deposi- 
tion from a jet, and is overpressured relative to 



the ambient ISM. Monitoring the intensity and 
searching for proper motions in the knots of the 
weak counterjet on this side of the source would 
provide important diagnostics for establishing the 
current state of the jet within the southwest lobe. 



• To order of magnitude, the mass in the hot shell is 
consistent with material swept up from the ISM. 
There must have been sufficient time since any 
earlier epoch of lobe expansion in the region for 
the ISM to be replenished. 



• The shell is overpressured compared with the 
minimum-energy pressure in the radio lobe (in 
magnetic field and radiating electrons) by a factor 
of ~ 10. There is no particular reason to think that 
a dynamical object, seen in a snapshot, should be 
in a state of minimum energy. However, if we do 
assume minimum energy in the lobe, and that the 
shell has reached equilibrium [but note that the 
sound-crossing time in the shell (thickness ~ 0.3 
kpc, c ~ 9 x 10~ 7 kpc yr -1 ) is about 15 per cent 
of the maximum time we estimate it has taken the 
lobe to reach its current size], the shell's overpres- 
sure relative to the radio lobe could be balanced 
by the ram pressure from internal motions in the 
lobe for a moderate relativistic proton loading. 



The shell's kinetic energy is ~ 5 times its ther- 
mal energy, and exceeds the thermal energy of the 
ISM within 15 kpc of the centre of the galaxy. 
As the shell dissipates, most of the kinetic energy 
should ultimately be converted into heat and this 
will have a major effect on Cen As ISM, provid- 
ing distributed heating. 
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